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Abstract 

In this paper, we study the entanglement properties of a spin-1 model 
the exact ground state of which is given by a Matrix Product state. The 
model exhibits a critical point transition at a parameter value a = 0. The 
longitudinal and transverse correlation lengths are known to diverge as a — ► 
0. We use three different entanglement measures S(i) (the one-site von 
Neumann entropy), S(i,j) (the two-body entanglement) and G(2,n) (the 
generalized global entanglement) to determine the entanglement content of 
the MP ground state as the parameter a is varied. The entanglement length, 
associated with S(i,j), is found to diverge in the vicinity of the quantum 
critical point a = 0. The first derivative of the entanglement measure E 
(= S(i), S(i,j)) w.r.t. the parameter a also diverges. The first derivative of 
G(2,n) w.r.t. a does not diverge as a — > but attains a maximum value at 
a = 0. At the QCP itself all the three entanglement measures become zero. 
We further show that multipartite correlations are involved in the QPT at 
a = 0. 

I. INTRODUCTION 

The entanglement characteristics of the ground states of many body Hamiltoni- 
ans describing condensed matter systems constitute an important area of study 
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in quantum information theory. Entanglement is an essential resource in quatum 
computation and communication protocols. Condensed matter, specially, spin sys- 
tems have been proposed as candidate systems for the realization of some of the 
protocols. Entanglement provides a measure of non-local quantum correlations in 
the system and it is of significant interest to determine how the correlations associ- 
ated with the ground state of the system change as one or more than one parameter 
of the system is changed. The focus on ground state characteristics arises from 
the possibility of quantum phase transitions (QPTs) which occur at temperature 
T = (when the system is in its ground state) and are driven solely by quantum 
fluctuations [T|. A QPT is brought about by tuning a parameter, either external 
or intrinsic to the Hamiltonian, to a special value termed the transition point. In 
thermodynamic critical phenomena, the thermal correlation length diverges and 
the thermodynamic quantities become singular as the critical point is approached. 
In the quantum case, the correlation length diverges in the vicinity of the QCP and 
the ground state properties develop non-analytic features. An issue of considerable 
interest is whether the quantum correlations, like the usual correlation functions, 
become long-ranged near the QCP. In a wider perspective, the major goal is to 
acquire a clear understanding of the variation in entanglement characteristics as a 
tuning parameter is changed. QPTs have been extensively studied in spin systems 
both theoretically and experimentally. In recent years, several theoretical studies 
have been undertaken to elucidate the relationship between QPTs and entangle- 
ment in spin systems [3], HJ El El [7]. In particular, a number of entanglement 
measures have been identified which develop special features close to the transition 
point. One such measure is concurrence which quantifies the entanglement between 
two spins (S = |). At a QCP, as illustrated by a class of exactly-solvable spin 
models (S = ~), the derivative of the ground state concurrence has a logarithmic 
singularity though the concurrence itself is non- vanishing upto only next-nearest- 
neighbour-distances between two spins [SIE]. Discontinuities in the ground state 
concurrence have been shown to characterize first order QPTs [HI M, ID]. Later, 
Wu et al. [5] showed that under some general assumptions a first order QPT, 
associated with a discontinuity in the first derivative of the ground state energy, 
gives rise to a discontinuity in a bipartite entanglement measure like concurrence 
and negativity. Similarly, a discontinuity or a divergence in the first derivative of 
the same entanglement measure is the signature of a second order phase transition 
with a discontinuity or a divergence in the second derivative of the ground state 
energy. Another measure of entanglement, studied in the context of QPTs, is the 
entropy of entanglement between a block of L adjacent spins in a chain with the 
rest of the system [4|. At the QCP, the entropy of entanglement diverges logarith- 
mically with the length of the block. There is, however, no direct relation with 
the long range correlations in the system. 
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A number of entanglement measures have recently been proposed which are 
characterized by a diverging length scale, the entanglement length, close to a QCP. 
The localizable entanglement (LE) between two spins is defined as the maximum 
average entanglement that can be localized between them by performing local 
measurements on the rest of the spins [TT]. The entanglement length sets the scale 
over which the LE decays. The two-body entanglement S(i,j) is a measure of the 
entanglement between two separated spins, at sites i and j, and the rest of the 
spins [7]. Let p(i,j) be the reduced density matrix for the two spins, obtained 
from the full density matrix by tracing out the spins other than the ones at sites i 
and j. The two body entanglement S(i,j) is given by the von Neumann entropy 

S(i,j) = -Trp(i,j)log 2 p(i,j) (1) 

In a translationally invariant system, S depends only on the distance n =| j — i \. 
As pointed out in [7], the spins that are entangled with one or both the spins at 
sites i and j contribute to S. The following results have been obtained in the 
case of the S = | exactly solvable anisotropic XY model in a transverse magnetic 
field. The model, away from the isotropic limit, belongs to the universality class 
of the transverse Ising model. The two-body entanglement S(i,j) has a simple 
dependence on the spin correlation functions in the large n limit. Away from the 
critical point, S(i,j) is found to saturate over a length scale as n increases. 
Near the QCP, one obtains 

S(i,j)- S(oo) ~n- 1 e"e^ (2) 

The entanglement length (EL), has an interpretation similar to that in the 
case of LE. The EL diverges with the same critical exponent as the correlation 
length at the QCP. S(i,j) thus captures the long range correlations associated 
with a QPT. At the critical point itself, S(i,j) — S(oo) has a power-law decay, i.e., 
S(i,j) — S(oo) ~ n~2. In the limit of large n, the first derivative of S(i,j) w.r.t. a 
Hamiltonian parameter develops a A— like cusp at the critical point. The univer- 
sality and a finite-size scaling of the entanglement have also been demonstrated. 
The one-site von Neumann entropy 

S(i) = -Tr p(i) log 2 p{i) (3) 

is also known to be a good indicator of a QPT [3j. It provides a measure of how a 
single spin at the site i is entangled with the rest of the system. The reduced density 
matrix p{i) is obtained from the full density matrix by tracing out all the spins 
except the one at the site i. Oliveira et al. [6] have proposed a generalized global 
entanglement (GGE) measure G(2, n) which quantifies multipartite entanglement 
(ME). G(2,n) for a translationally symmetric system is given by 



3 



G ( 2 ' n ) = J^T [1 " S I l>0',J +")]im I 2 ] (4) 

l,m=l 

where p(j,j + n) is the reduced density matrix of dimension d. The factor 2 in 
(7(2, n) indicates that the reduced density matrix is that for a pair of particles. Wu 
et al. [5] considered QPTs characterized by non-analyticities in the derivatives of 
the ground state energy. These arise from the non-analyticities in one or more of 
the elements of the reduced density matrix. In terms of the GGE, a discontinuity 
in (7(2, n) signals a first order QPT, brought about by a discontinuity in one or 
more of the elements, [pj,j+ n ]lm of the reduced density matrix [6j. A discontinuity 
or divergence in the first derivative of G(2,n) w.r.t. the tuning parameter occurs 
due to a discontinuity or divergence in the first derivetives of one or more of the 
elements of the reduced density matrix. The associated QPT is of second order. 
Non-analyticities in (7(2, n) and its derivatives thus serve as indicators of QPTs. 
In the case of the XY S — | spin chain, the GGE measure shows a diverging EL 
as the QCP is approached. The EL = ^ where £c* is the usual correlation 
length. Thus, both the length scales diverge with the same critical exponent near 
the QCP. 

The relationship between entanglement and QPTs has mostly been explored for 
spin-| systems. The entanglement properties of the ground states of certain spin— 1 
Hamiltonians have been studied using different measures [HI H21 [13] . Numerical 
studies show that the LE has the maximal value for the ground state of the spin-1 
Heisenberg antiferromagnet with open boundary conditions (OBC) [33]. In the 
case of the spin-1 Affleck-Kennedy-Lieb-Tasaki (AKLT) model [H], the result can 
be proved exactly. A class of spin-1 models, the 0-deformed AKLT models, is 
characterized by an exponentially decaying LE with a finite EL The length 
£e diverges at the point = though the conventional correlation length remains 
finite [13]. A recent study [15] shows that in the case of spin-1 systems, the use 
of LE for the detection of QPTs is not feasible. An example is given by the S = 1 
XXZ Heisenberg antiferromagnet with single-ion anisotropy. The model has a rich 
phase diagram with six different phases. The LE is found to be always 1 in the 
entire parameter region and hence is insensitive to QPTs. The ground states of 
certain spin-1 models have an exact representation in terms of matrix product 
states (MPS) [HI [TZI [TH]. The ground state of the spin-1 AKLT model, termed 
a valence bond solid (VBS) state, is an example of an MPS. The ground state 
is characterized by short-ranged spin-spin correlations and a hidden topological 
order known as the string order. The excitation spectrum of the model is further 
gapped. In the MPS formalism, ground state expectation values like the correlation 
functions are easy to calculate. This has made it particularly convenient to study 
phase transitions in spin models with MP states as exact ground states [T7] . The 
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transitions identified so far include both first and second order transitions and 
are brought about by the tuning of the Hamiltonian parameters. The second 
order transition in the class of finitely correlated MP states, however, differs from 
the conventional QPT in one important respect. The spin correlation function 

_ n 

is always of the form A c e ? c for large n. The correlation length £c* diverges as 
the transition point is approached. The pre-factor Ac, however, vanishes at the 
transition point [17]. This is in contrast to the power-law decay of the correlation 
function at a conventional QCP. Some distinct features of QPTs in MP states have 
recently been identified [19J. One of these relates to the analyticity of the ground 
state energy density for all values of the tuning parameter. In a conventional QPT, 
the energy density becomes non-analytic at the QCP. The MP states appear to 
provide an ideal playground for exploring novel types of QPTs. In this paper, 
we consider a spin-1 model, the exact ground state of which is given by an MP 
state [20]. The model has a rich phase diagram with a number of first order phase 
transitions and a critical point transition. We study the entanglement properties 
of the ground state with a view to pinpoint the special features which appear 
close to the critical point. This is done by using three different entanglement 
measures, namely, the single-site, two-body and generalized global entanglement 
defined earlier. 

II. REDUCED DENSITY MATRIX OF MP 
GROUND STATE 

We consider a spin-1 chain Hamiltonian proposed by Klumper et al. [20] which 
describes a large class of antiferromagnetic (AFM) spin-1 chains with MP states 
as exact ground states. The Hamiltonian satisfies the symmetries : (i) rotational 
invariance in the x — y plane, (ii) invariance under S z — > —S z and (iii) translation 
and parity invariance. The Hamiltonian has the general form 

L 
3=1 

h j>j+1 = a A 2 j + a 1 {A j B j + BjAj) + a 2 B) + a 3 Aj + a 4 Bj(l + Bj)+ 

+ a 5 ((Stf + (S* +1 f + C (5) 

where L is the number of sites in the chain and periodic boundary conditions (PBC) 
hold true. The parameters aj are real and C is a constant. The nearest-neighbour 
(n.n.) interactions are 
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a . — qx qx _i_ qy qy 

The constant C in Eq. (5) may be adjusted so that the ground state eigenvalue of 
hjj + i = 0. Hence 

h jtj+1 >0 ^H>0 (7) 

i.e., H has only non-negative eigenvalues. In the AFM case, the z-component of 
the total spin of the ground state S^ ot = 0. KKimper et al. showed that in a 
certain subspace of the <x, —parameter space , the AFM ground state has the MP 
form. Let |0) and |±) be the eigenstates of S z with eigenvalues 0, +1 and — 1 
respectively. Define a 2 x 2 matrix at each site j by 

with non- vanishing parameters a, a ^ 0. 
The global AFM state is written as 



\tpo (a, a)) = Tr ( gi <g> g 2 <g> <g> g L ) (9) 

where '®' denotes a tensor product. One can easily check that S* ot j^o) = 0, i.e., 
the state is AFM. One now demands that the state \if) (a, a)) is the exact ground 
state of the Hamiltonian H with eigenvalues 0. For this, it is sufficient to show 
that 

h 3,i+i (9j ®9j+i) = (10) 
Eq. (3) and (10) are satisfied provided the following equalities 

1) cr = sign(a 3 ), 2)aa = a 3 -a 1 , 

3) «5 =| «3 | +Q!o(l — a 2 ), 4) a.2 = a^a 2 — 2 | a \ 

and inequalities 

a ^ 0, a 3 7^ 0, a 4 > 0, a > (12) 

hold true. The state j^o (a, a)) is the ground state of the Hamiltonian (5) with 
ground state energy zero provided the equalities in (8) are satisfied. The inequal- 
ities constrain the other eigenvalues of hjj + i to be positive. If the inequalities 
are satisfied, the ground state can be shown to be unique for any chain length 
L. Also, in the thermodynamic limit L — > oo, the excitation spectrum has a gap 
A. With equality signs in the inequalities (12), the state \ipo (a, a)) is still the 
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ground state but is no longer unique. The spin— 1 model has the typical feature 
of a Haldane-gap (HG) antiferromagnet. In fact, the AKLT model is recovered as 
a special case with a — 2, & = 1, «3 = Sao > 0, a 2 = — 2a and = 3«o- The 
state (9) now represents the VBS state. 

Using the transfer matrix method [16], the ground state correlation functions 
can be calculated in a straightforward manner. The results are (L — > oo, r > 2) : 
Longitudinal correlation function 

< s ^> = -(T^(ttr)' (13) 

Transverse correlation function 

(S*S?) = -\a\ [a + sign a] ( 14 ) 

The correlations (13) and (14) decay exponentially with the longitudinal and trans- 
verse correlation lengths given by 

er^Mi + H) (15) 

Furthermore, the string order parameter has a non-zero expectation value in the 
ground state. One finds that the correlation lengths diverge as a — > 0. At the 
point a = 0, the correlation functions given by Eq. (13) and (14) are zero. At 
a conventional QCP, the correlation functions have a power-law decay. We will, 
however, refer to the point as a QCP since the correlation lengths diverge as the 
point is approached. A consequence of the diverging correlation length is that the 
excitation spectrum of the spin-1 model, which is gapped (the Haldane phase) for 
a > 0, becomes gapless at the critical point a = [20]. The presence or absence of 
a gap in the excitation spectrum of a system is reflected in the low temperature 
thermodynamic properties of the system. Furthermore, the string order parameter 
has a non-zero expectation value in the ground state for a > and becomes zero 
at a = indicating the appearance of a new phase. Refs. [TH [T7] provide several 
other examples of spin-1 models with finitely correlated MP states as exact ground 
states. All these models exhibit critical point transitions with features similar to 
those in the case of the spin-1 model described by the Hamiltonian in Eq. (5). 
We now focus on the entanglement properties of the MP ground state (Eq. (9)). 
We consider a to be > and a = +1 in Eq. (8). The one-site reduced density 
matrix p(i) (Eq. (3)) obtained by tracing out all the spins except the i-th spin 
from the ground state density matrix p = |^ )(' i / ; o|, can be calculated using the 
transfer matrix method [TB]. The density matrix, from Eq. (9), is 

(16) 

{n a ,m a } 



—1 



In 



1 + lal 



1 — a 
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The summation is over all the indices, rij, m«, % = 1,2, L. 

We define a 4 x 4 matrix / (the elements of which are operators) at any lattice 
site as 

//il/U2 =^ /(ni,mi)(ri2,m2) = 9n\ni drrixm? (l^) 

The convention of the ordering of the multi-indices is // = 1, 2, 3, 4 <-> (11), (12), (21), (22). 
Thus, / can be written as 



/ |0> <0| 



/ 



a|0) (1| 



a |1> <0j a|l)(l| \ 



v / 2|0)(-l| — 10> <0| -a|l)(-l| v^|l) 



a | — 1> <0| —a |— 1> <1| 



-|o> (0| 



a|0) (1| 



V a|— 1> <— 1| -v^|-l)(0| -v^|0)(-l| |0)(0| / 



Also, 



(18) 



(19) 



pW = ^UI^o) (Vol 

where the trace is over all the spins except the i-th one. The transfer matrix F at 
a site m is obtained by taking the trace over / at the same site, i.e., 



F m = J2 ( k \ fm \k) 



(20) 



where the states \k) are the states |0), |±1). The transfer matrix F is obtained as 

(21) 



/ 1 a \ 
0-100 

0-10 

V a 1/ 



The eigenvalues are 

ex = 1 + a, e 2 = I - a, e 3 = -1, e 4 = -1 
The corresponding eigenvectors are 



(22) 



e 3 



/ 1 \ 




V i J 
1 



V o / 



1^ = 73 



f ~ 1 \ 




V i / 



1 

V o / 



(23) 



From Eq. (20), 



J2a=l ( e a 


F L-lf 






a 


F L 





(24) 



The factor in the denominator takes care of the condition Tr p = 1. On taking the 
thermodynamic limit L — > oo, we get 



p(i) = £ / (ei|/|ei) 
In the |0, ±1) basis, the reduced density matrix becomes 



(25) 



/ -p— \ 

l+a 

2(l+a) 
2(Tfa) ) 



\ 



(26) 



The calculation of the two-site reduced density matrix p(i,j) follows in the 
same manner. p(i,j) is given by 



p(ij) = Tr{ 3 L |V>o) <-0oj 
where the trace is taken over all the spins except the i-th and j-th ones. 

P(hj) = 

In the thermodynamic limit L — > oo, p(i,j) reduces to 

p(i,j) = Y, £ » 2 (-r) ( e il / l e «) ( e °l / l e i) 

a=l V£ l y 





f pj-i- 


1 f pL-j 


e a ) 









(27) 



(28) 



(29) 



where n = \j — i\. 

The matrix p(i,j) is a 9 x 9 matrix and defined in the two-spin basis states 
\lm) with the ordering 

\lm) = |11) , 1 10) , |01) , |1 - 1) , |-11) , |00) , |0 - 1) , |-10) , |-1 - 1) (30) 

The non-zero matrix elements, b pq (p = 1, ...,9, q — 1, ...,9), of p(i,j) are : 

h -U _ _a? « 2 fi-a\ n 

"11 - "99 - 4(i +a )2 4(l-a 2 ) \l+a J 

^22 = &33 = ^77 = &88 = 2(l+a) 2 

b 4A = b 55 = 4(1 + a)2 + 4(1 ° o2) (i^f) (31) 

&23 = ^32 = ^46 = &64 = &56 = ^65 = ^78 = ^87 = 2 (l+a) (~ l+o J 



^66 - (T+^)2 

It is easy to check that p(i,j) has a block-diagonal form. 
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III. ENTANGLEMENT MEASURES S(i), S{iJ), 
G(2,n) 

We now determine the entanglement content of the ground state \ipo) (Eq. (9)) 
of the Hamiltonian (Eq. (5)) using the entanglement measures S(i), S(i,j), and 
G(2,n). The calculations are carried out for different values of the parameter a 
in Eq. (8). The ultimate aim is to probe the special features, if any, of entangle- 
ment in the vicinity of the QCP at a = 0. From Eq. (3) and (26), the one-site 
entanglement 

S(i) = Y^K 1 + a ) l °92 (1 + a) - a log 2 a + a] (32) 

Figure 1 (top) shows the variation of S(i) w.r.t. a. The one-site entanglement has 
the maximum possible value log 2 3. This is attained at the AKLT point a = 2. The 
VBS state is in this case the exact ground state. In the VBS state, each spin-1 at 
a specific lattice site can be considered as a symmetric combination of two spina's 
[U]. In the VBS state, each spin-| at a particular lattice site forms a spin singlet 
with a spin-i at a neighbouring lattice site. S(i) has the value zero at the QCP 
a = 0. Figure 1 (bottom) shows the variation of ^§j^ with the parameter a. The 
derivative diverges as the QCP is approached. This is the expected behaviour at 
the QCP of a conventional QPT. In the latter case, however, S(i) has the maximum 
value at the QCP [3]. 

From Eq. (1) and (31), the two-body entanglement S(i,j) is 

9 

S(i,j) = -J2\ilog 2 \i (33) 

Where Aj's are the eigenvalues of the reduced density matrix p(i,j). These are 
given by 

\ _ \ _ a a ( ]_\ n 

Ai — A 2 — 2 (l+o)2 2(l+o) V 1+aJ 

A 3 = A 4 = 2(1 + a)2 + 2 (l+a) 

\ a 2 I a 2 ( 1— a \ n 

A 5 — 4(l +0 )2 -T 4(i_ a 2) \l +a ) 

x _ x _ __g? af_ (l-a\ n 

Aq— A? — 4(l+o) 2 4 (l-a 2 ) U+a/ 

^ 8 = 2 (4(l+a) 2 + 4(l-a 2 ) (l+f) + TT^ 2 ) ~~ 2(l+a) (^4) 

(a 2 ~4) 2 , n n 2 ( l_\ 2n , a 4 / l^aA 2 " , a 2 (a 2 -4) /l_a\ J 

16(l+a) 2 "r Z " V 16(l-a) 2 V 1+"/ 8(l-a 2 ) V 1+a / 



\ _ 1 ( a 2 , a 2 / l-a\" , 1_ 

A 9 — 2 V4(l+a) 2 4(l-a 2 ) V 1+a/ 1+a 



1 



1+a 2 J ^ 2(l+a) 

f (a 2 -4) 2 9 2 / ^\ 2 « a 4 /l^ 2 ™ , a 2 (a 2 -4) /W ^ * 

I 16(l+a) 2 ^ ZU V 1+a/ ^ 16(l-a) 2 V 1+a/ 8(l-a 2 ) V 1+0 
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Knowing the reduced density matrix p(i,j), the correlation functions (S" S"} 
(a = x, y, z) can be calculated in the usual manner. One then recovers the ex- 
pressions in Eq. (13) and (14) (r = n + 1, where n = \j — Figure 2 (top) 
shows the variation of S(i,j) as a function of a for n = 1000. Figure 2 (bottom) 
shows the variation of the derivative ^ w.r.t. a for the same value of n. The 

aa 

maximum of S(i,j) is at the AKLT point a = 2 and has the value zero at a = 0. 
For large n, the derivative 9S g a ^ diverges near the QCP at a = 0. The last fearure 
is characteristic of a conventional QPT [7]. 

We next calculate the GGE G(2, n) (Eq. (4)). This is easily done as the matrix 
elements of the reduced density matrix (Eq. (3)) are known. Figure 3 (top) shows 
the variation of G(2, n) versus a for n = 1000. Figure 3 (bottom) shows the plot of 
9G ^ n ^ against a. Again G(2, n) has the maximum value at the AKLT point and is 
zero at a = 0. The derivative 9G ^ n ^ does not diverge as a — > in contrast to the 
case of a conventional QPT [6j. The derivative, however, attains the maximum 
value at the QCP a = 0. Figure 4 (top) shows the plots of S(i), S(i,j), and G(2, n) 
versus a for n = 1000. Figure 4 (bottom) shows the plots of the first derivatives 
of the same quantities w.r.t. a for n = 1000. The plots are shown for comparing 
the different entanglement measures. 

We next determine the EL and its variation w.r.t. the parameter a. We 
consider the entanglement measure S(i,j) for this purpose. Close to the QCP 
a = and in the limit of large n, one can write 

_ n 

S(n = \j-i\)- S{oo) ~ A e e *S (35) 



The longitudinal and transverse correlation functions, p z n = (^SlS^ +l j and p* = 

(SfS* + i} are given by Eq. (13) and (14) with r = n + 1. For a < 1, p z n decays 
faster than p x n with n. The eigenvalues Aj's, i = 1, ...9, can be expressed in terms 
of the correlation functions p z n and p x n . For large n, the contributions from p z n can 
be ignored. The eigenvalues now reduce to the expressions 

Ai = A 2 = 2 {i+ a y ~ ^Pn 

A 3 = A 4 = 2 (l+a)2 + 4 Pn 

\ — °- 2 
A 5 — 4 (i +a )2 

A fi = Ay = q2 , (36) 



4(l+o) 

\ _ 1 ( a 2 , _J_A _ 1 / (° 2 ~4) 2 _j_ o („ x\2^ 
A § — 2 V4(l+a) 2 1+a 2 J 2 V 16(l+a) 2 ^ Z 



n 

\ - 1 ( a2 I M I if ( q2 ~ 4 ) 2 I Q rn^^ 2 ^ ^ 
A 9 — 2 V4(l+a) 2 1+a 2 / 2 V 16(l+a) 2 ^ Z W 

From Eq. (1) and (36), one ultimately arrives at the expressions 



S(n = \j-i\)- S(oo) ~ A e (p x n ) 2 ~ A e e"^ (37) 
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The pre-factor A e = at the QCP a = . The EL £g is given by 

& = 2 = 2Zn(l + a) (38) 

where £ t is the transverse correlation length (Eq. (15)). In the case of the S = | 
anisotropic XY model in a transverse field, an expression similar to that in Eq. 
(37) is obtained close to the QCP in the limit of large n [7]. The pre-factor in this 
case, however, does not vanish at the QCP but has a power-law dependence on n. 
Figure 5 shows the variation of £e w.r.t. a based on the entanglement measure 
S(i,j). 

The total correlations, with both classical and quantum components, between 
two sites i and j are quantified in terms of the quantum mutual information [2T1 [22] 



h^Sty + SW-SiiJ) (39) 

As explained in [21], a comparison of the singular behaviour of S(i) with that of 
Iij allows one to determine whether two-point (Q2) or multipartite (QS) quantum 
correlations are important in a QPT. Figure 6 (top) shows a plot of versus a 
for n = 1000. Figure 6 (bottom) shows the variation of versus a for n = 1000. 
The derivative does not diverge as a — > 0, a behaviour distinct from that of 
close to a = 0. The difference in the singular behaviour of quantities associated 
with S(i) and iy indicates that multipartite quantum correlations are involved in 
the QPT. 



IV. DISCUSSIONS 

The MP states provide exact representations of the ground states of several spin 
models in low dimensions [32]. The remarkable features of such states arises from 
the fact that complicated many body states have a simple factorized form. The 
simplicity in structure makes the calculation of the ground state expectation values 
particularly easy to perform. The spin-1 AKLT model is a well-known example 
of a spin model in Id the exact ground state of which (a VBS state) has an MP 
representation. The AKLT model and the spin-1 Heisenberg AFM chain belong 
to the same universality class [23]. The insight gained from the study of models in 
the MP formalism is expected to be of relevance in understanding the properties 
of more physical systems. The MP states also serve as good trial wave functions 
for standard spin models. The MP representation lies at the heart of the powerful 
density matrix renormalization group (DMRG) method and provides the basis for 
several interesting developments in quantum information [24]. 

Studies of the entanglement characteristics of the MP states have begun only 
recently. The QPTs which occur in such states have characteristics different from 
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FIG. 3: Plot of G(2, n) (top) and (bottom) as a function of a. 
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FIG. 4: Plots of S(i) (q), S(i,j) (r), and G(2,n) (p) (top) and the corresponding 
first derivatives w.r.t. a (bottom) as a function of a. E represents the entanglement 
measure. 
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FIG. 5: Plot of EL as a function of a. 
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FIG. 6: Plots of hj (top) and (bottom) as a function of a. 
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those of conventional QPTs. It is thus of considerable interest to determine whether 
the entanglement content of MP states develops special features close to a QCP. 
In this paper, we consider a spin-1 model the exact ground state of which is of the 
MP form over a wide range of parameter values. The model exhibits a novel QPT 
in that the longitudinal and transverse correlation lengths diverge as the QCP is 
approached but the correlation functions vanish at the QCP. In a conventional 
QPT, the correlation functions have a power-law decay at the QCP. In the spin-1 
model, the divergence of correlation lengths is accompanied by the excitation gap 
going to zero. The string order parameter, which has a non-zero expectation value 
in the MP state for a > 0, vanishes at the QCP a = 0. The distinctive signatures 
indicate the appearance of a new phase. We study the entanglement properties 
of the MP state for different values of the parameter a. The measures used are 
S(i) (one-site von Neumann entropy), S(i,j) (two-body entanglement) and G(2, n) 
(GGE). All the entanglement measures have zero value at the QCP so that the 
ground state is disentangled at that point. As seen from the different plots, figures 
(1), (2), (3), and (4), the entanglement content, as measured by E =S{i), S(i,j) 
and G(2,n), has a slow variation w.r.t. a for a > 2. At the AKLT point a = 2, 
E reaches its maximum value and as a is reduced further, the magnitude of E 
falls rapidly to approach zero value at a = 0. The study of conventional QPTs 
shows that E is maximum at a QCP [H EE [7]. Also, ^ diverges as the QCP is 
approached. The EL, as calculated from S(i,j) and G(2,n) for large n, also 
diverges with = & where £c is the usual correlation length. In the case of 
the spin-1 model under consideration, ^ diverges as a — > when E = S(i) and 
S(i,j). The EL, as calculated from S(i,j) in the large n limit, also diverges 
with £e = ^2~. One now has the interesting situation that the entanglement content 
of the MP state decreases as a — > but the entanglement is spread over larger 
distances. The derivative 9G ^ n ^ , however, does not diverge as a — > but attains 
a maximum value at the QCP. The results can be understood by noting that in 
all the three cases, E = S(i), S(i,j) and G(2,n), the reduced density matrices 
p(i) and p(i,j) smoothly approach the forms associated with pure states as the 
parameter a tends to zero. The matrix elements of the reduced density matrices do 
not develop non-analyticities in the parameter region of interest. Thus, the energy 
density, calculated from the reduced density matrix p(i,j), does not develop a non- 
analyticity at the QCP. The derivative aG ^'"- > depends upon the first derivatives 
of the matrix elements of p(i,j). Since the latter is analytic for a > 0, 8G ^' n - > does 
not diverge in the whole parameter regime including the point a = 0. In the cases 
of the entanglement measures E = S(i) and S(i,j), the derivative diverges as 
a — > due to the divergence of log 2 a in the same limit. The divergence is thus due 
to the special form of the von Neumann entropy and occurs for n > 1. A recent 
work [25] provides another example of such a singularity. Though dG Q^ does 
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not diverge or become discontinuous at the QCP a = 0, it attains its maximum 
value at the point. The first derivative of the string order parameter w.r.t. a also 
attains its maximum value at a = though the order parameter itself vanishes 
at the point. Figure (4) shows that amongst the three entanglement measures 
E =S(i), S(i,j) and G(2, n), used in this study to obtain a quantitative estimate 
of the entanglement content of the MP ground state, the measure S(i,j) yields 
the largest value of the entanglement at different values of a. The difference in 
the singular behaviour of the measures S(i) and the mutual information entropy 
7jj as a — > indicates that multipartite quantum correlations are involved in the 
QPT. In summary, the present study identifies the entanglement characteristics 
of the MP ground state of a spin-1 model close to the critical point a = 0. The 
features are distinct from those associated with conventional QPTs. Several spin 
models are known for which the MP states are the exact ground states [T7]. Some 
of these models have interesting phase diagrams exhibiting both first and second 
order phase transitions. It will be of interest to extend the present study to other 
spin models (both S = | and 1) in order to identify the universal characteristics 
of QPTs in MP states. 
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